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I. INTRODUCTION 

In recent years, the properties of the inner crust of neu- 
trons stars have been investigated vifithin various frame- 
virorks, especially focusing on their microscopic struc- 
ture and superfluid properties. The microscopic calcu- 
lations of the inner crust matter are usually studied in 
the Wigner-Seitz (WS) approximation [Ij. Following 
the standard approaches, the inner crust consists of a 
lattice of WS cells, each cell containing a neutron-rich 
nucleus immersed in a sea of dilute gas of neutrons and 
relativistic electrons uniformly distributed inside the cell 
Q.The first microscopic calculations of the properties of 
the inner crust matter were done by Negele and Vau- 
therin in the 70s The calculations were performed 
assuming a set of non-interacting cells described in the 
Hartree-Fock (HF) approach with 11 representative cells. 
These cells are distributed in different zones of the inner 
crust with densities covering a range from 1.743 x 10"'^ po 
to 0.5yOo, Po = 0.16 fm~^ being the nuclear matter satu- 
ration density. The optimal number (N, Z) of neutrons 
and protons in each cell is obtained by searching for the 
lowest binding energy satisfying the /3-stability condition 
of the cell. 

More recently, the superfluid properties and their in- 
fluence on the specific heat were investigated in the 
self-consistent Hartree-Fock-Bogoliubov (HFB) approach 
[1, The collective excitations and the cooling time 
of the inner crust of neutron stars were also studied in 
the framework of a spherically symmetric HFB + quasi- 
particle random phase approximation [s], Q and of HFB 
approach at finite temperature Q- 

All the above calculations were done within the non- 
relativistic framework with the SLy4 Q Skyrme interac- 
tion in the mean field channel and a density-dependent 
delta force for the pairing interaction. Later on, the 
properties of the WS cells have been studied in the non- 
relativistic work of Baldo et al. [ol-flU using an energy 
functional method including the pairing correlations of 
protons and neutrons. The relativistic Hartree, or Rel- 
ativistic Mean Field (RMF) approach [l2| has also been 
applied to investigate the structure of WS cells and the 



influence of different boundary conditions on the struc- 
ture of the cells. 

However, there is so far no study of the inner crust 
in the HF approach using flnite range interactions. The 
aim of this paper is to examine the behaviour of density 
distributions of protons and neutrons in the different WS 
cells obtained with a flnite range interaction within a self- 
consistent HF-BCS calculation, where the HF equations 
are solved with the Dirichlet-Neumann mixed boundary 
conditions. For the numerical studies one obvious choice 
is the DIS interaction [isj which is widely employed in 
flnite nuclei calculations. More recently, the M3Y-P4 in- 
teraction proposed by Nakada [l3| seems to work satis- 
factorily well for describing nuclear ground states in the 
HFB framework. This flnite-range effective interaction 
derives from the well-known M3Y interaction and it is 
interesting to study its predictions for neutron stars. 

The structure of this paper is as follows. In Section II 
we present the HF and HF-BCS approaches using flnite- 
range density-dependent interactions in both the mean 
fleld and pairing channels. The results of the calcula- 
tions arc presented in Section HI. The evolution of the 
neutron and proton densities as well as the neutron pair- 
ing flclds in the regions of the inner crust are discussed. 
The summary and perspectives of this study are given in 
Section IV. 



II. HF-BCS CALCULATIONS WITH FINITE 
RANGE INTERACTIONS 

A practical difficulty of HF calculations in WS cells 
is that the cell radius can be large (up to 40-50 fm). 
Solving directly the integro-differential equations in co- 
ordinate space can lead to inaccuracies and instabilities 
[l5| . The method that we adopt consists in solving the 
HF equations in a basis well adapted to each WS cell un- 
der study. In subsection A we derive the HF equations 
in coordinate space, and in subsection B we give some 
details on their solution. 
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A. HF equations in coordinate space 



where 



We summarize here the analytical expressions needed 
for HF calculations with a finite-range interaction, ta king 
as an illustrative example the Gogny interaction [l3l Il6| 
which contains a sum of two Gaussians, a zero-range 
density-dependent part and a zero-range spin-orbit part. 
Expressions for other types of finite-range interactions 
can be easily deduced from this case. For example, 
the M3Y-Pn interactions of Nakada have zero-range, 
density-dependent terms just like Gogny and Skyrme in- 
teractions and therefore, the corresponding contributions 
to the mean field potentials can be deduced from the ex- 
pressions given here. In addition, the M3Y-Pn interac- 
tions also contain tensor terms, except for the M3Y-P4 
interaction that we use in this work and therefore, we do 
not discuss the tensor contributions. Finally, the original 
M3Y-P4 interaction contains short-range spin-orbit com- 
ponents (with fi = 0.25/m and 0.40/to, see Eq. [T]) that 
we approximate for simplicity by a zero-range spin-orbit 
component as in the Gogny force. The corresponding 
spin-orbit strength Wq = IQOAleV.fm^ is adjusted on 
the empirical lpl/2 — lp3/2 proton splitting in ^^O as in 
[l6| . Throughout this work we limit ourselves to spheri- 
cally symmetric systems, which should apply to baryonic 
densities in the inner crust ranging from 1.410~'^po to 
about 0.5/90, where pQ=0.16/m~'^ is the nuclear matter 
saturation density [13 . 

In order to obtain the Hartree and Fock potentials, we 
use the multipole expansion of the Gaussian form factor 
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We express the central part of the Gogny force in the 
form 



SLJ v=l 
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A,(5 = 0) = W,-H,P^^ 
B — M P"^ 



(3) 



Here, W^, B^, Hi,, M^, are the parameters of the 
Gogny interaction or appropriate combinations of the pa- 
rameters of the M3Y-Pn interactions, P'^ and P"^ are the 
spin and isospin exchange operators, respectively. We 

have introduced the tensors y^"^^)*^ which are tensorial 

(m) 

products of a spherical harmonic Y^^ with a Pauli spin 
matrix 



(2) 



<SL)J ^^^s^Yt]. ^ (4) 

Because of the spherical symmetry assumption, the 
single-particle wave functions (fi{r, cr, q) can be factorizcd 
into a radial part Ui(r), a spin-angular part ^ijmi''^, u), 
and an isospin part (''')■ 

V'i{r,a,q) = -^^yij,n{f,(j)Xq{'r) (5) 

where 

■mi77ls 

(6) 

Xm, (c) being a spinor corresponding to a spin projec- 
tion uis, and the index i stands for the following set of 
quantum numbers: the charge q, the principal quantum 
number n, the orbital angular momentum Z, the total an- 
gular momentum j, and the magnetic quantum number 

TO. 

Using the tensors T^^j^^^^"^ and the single-particle wave 
hmctions ipi{r, cr, q) we can calculate the direct (Hartree) 
and exchange (Fock) potentials in each {I, j) partial wave. 
With the help of the results ([XB IM)) of Appendix A it 
is straightforward to obtain the radial integro-diffcrential 
HF equations in coordinate space for the radial wave 
functions Ui{ri): 



kik + i) 



+ + 1) - hik + 1) - t]M^, "('^i)"^(ri) = e^u,{n) 



(7) 



where the local central potential Uf contains the di- 
rect contributions of the finite range forces (nuclear and 
Coulomb) as well as the direct -I- exchange contributions 
of the zero-range density-dependent forces, while the di- 
rect -|- exchange contributions of the zero-range spin-orbit 



force are in both the Wg^^ component of U['{r) and in 
the one-body spin-orbit potential W^'^^Q: 

L/f (n) = C/f (n) + U^^ir,) + y^^(ri) + W^^\rr) . 

(8) 
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The non-local potential Uf is composed of the exchange 
contributions of all finite range (nuclear and Coulomb) 
forces: 

U^{ri,r2)^Ufin,r2) + Vf^iri,r2) . (9) 
The corresponding expressions are given in Appendix A. 

B. Solving HF equations in a basis representation 

Solving the HF equations directly in coordinate space 
has several advantages: the results do not depend on the 
choice of a basis and its truncation, and the individual 
wave functions have a correct asymptotic behavior. How- 
ever, in the regions of the inner crust where the density 
is below 10~^po the radius of the WS cells can be larger 
than 40 fm and strong numerical instabilities may appear 
for single-particle wave functions with large orbital mo- 
mentum (/ ~ 15 — 20). The cause of the problem is that 
the HF wave functions behave like r'+^ near the origin. 
Calculations requiring a certain number of nodes can fail 
if the HF potential is not accurate enough at a given iter- 
ation. To avoid this dificulty we have developed a basis 
expansion method [l^ for solving the HF equations as 
a matrix diagonalization problem. The advantage of the 
basis expansion method is that it produces a high ac- 
curacy (measured by the orthogonality of the solutions) 
with a smaller number of points, while this may not be 
obtained with the coordinate space method. 

The harmonic-oscillator functions basis is the most 
common choice for expanding the HF single-particle wave 
functions. Here, we choose instead, as a basis, the spher- 
ical Bessel functions normalized inside the Wigner-Seitz 
cell sphere. 

The radial part ""'^^''-^ of the single-particle wave func- 
tion in Eq. ([5]) can be expanded on the normalized spher- 
ical Bessel functions as 

^ = Ea,,j.(fcfV) (10) 

i—l 

where N is the dimension of the basis and ii{kf\) the 
normalized spherical Bessel functions with boundary con- 
ditions that will be specified later. The single-particle 
wave function </?i(r, cr, q) of Eq. ([5]) can be expanded as 

TV 

Vnljin{r,(y,q) ='^Cnlj,i'4^nl]m,i{'^,0,q) , (11) 
i=l 

where the orthonormal basis is 

V'«;jm,i(T',cr, (?) = j;(A:fV)y/j™(f,cr)xq(r) . (12) 

Under the / and j conservation the HF Hamiltonian ma- 
trix is decomposed into blocks: 

H'ij^ ^< ^llm.^\H\^lJm,^' > , (13) 



where H stands for the HF Hamiltonian. The HF equa- 
tions then become: 

N 
i' = l 

(11) 

where H-^, is a symmetric N x N matrix, 

HllP =|^(fcf^)''5,.+ (15) 
/ J Vif^in,r2)jiikl''>n)ji{kl!K2)dr^dr2 . 

Here, the potentials V^^^(ri,r2) can be obtained from 
the formulas (|A1I - IA6|) of Appendix A. The HF equations 
are solved by an iterative procedure, starting from an 
initial Woods-Saxon potential. In the present calculation 
we use a basis size N=20. The convergence criteria is set 
on the single-particle energies Cnijm to be 1 keV. 

The HF calculations are done with finite-range inter- 
actions (the DIS [ill and M3Y-P4 [Hi forces) imposing 
Dirichlet-Neumann boundary conditions at the edge of 
the cell as introduced in Ref. [i|. These boundary con- 
ditions for the single-particle wave functions are taken 
as follows: (i) the even parity wave functions vanish at 
the edge r = Rws of the box; (ii) the first derivative of 
the odd-parity wave functions vanish at r = Rws- The 
purpose of these chosen boundary conditions is to ob- 
tain an approximately constant density at large distance 
from the center of the cell, thus simulating a lattice of 
nucleus-like systems embedded in a uniform neutron gas. 

Of course, these boundary conditions are somewhat ar- 
bitrary and one could as well choose alternative boundary 
conditions. Thus, another kind of boundary conditions 
could be chosen in the following way: odd I wave func- 
tions vanish at r = Rws, the first derivatives of even I 
wave functions vanish at r = Rws- As shown in Ref. 
[m , the two kinds of boundary conditions can be used 
in the calculations of the neutron star inner crust. The 
difference of the binding energies per nucleon for each cell 
will increase with the increasing density (see Table 1 of 
Ref. [IH). However, the values of these uncertainties are 
smaller than the variations of the equilibrium configura- 
tion connected with the pairing effects [§, [ii| . Therefore, 
in our study the boundary conditions and the spherical 
WS cell structure are kept the same as in Ref. [l| . We do 
not redetermine the number (N, Z) of neutrons and pro- 
tons in the considered cells. This study can be envisaged 
in the near future. 

It should be noted that the Bloch boundary conditions 
can be used at the cost of more complex calculations as 
presented in Refs. [l"9l - [2l| . The validity in neutron star 
crust of the Wigner-Seitz approximation that we use here 
has been discussed in llil. 
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C. Pairing correlations in BCS approximation 

The effects of pairing correlations are known to be sub- 
stantial in the inner crust 0, Q and it is necessary to 
include them in order to obtain a more realistic picture. 
In this first study of the inner crust with finite range 
interactions in the particle-hole (mean field) channel we 
adopt the simplified BCS picture for describing the neu- 
tron pairing fields, but the HF-BCS model remains com- 
pletely self-consistent. Indeed, the BCS occupation num- 
bers are recalculated at each HF iteration. 

Various effective interactions can be used in the pairing 
channel. In this work the pairing correlations are treated 
in the BCS approximation with both the zero-range 
density-dependent and finite range density-dependent 
forces in the pairing channel. 

The importance of the density dependence of the pair- 
ing interaction is well known in the theories of superflu- 
idity in neutron stars. As shown in Ref. it is im- 
possible to deduce the magnitude of the pairing gaps in 
neutron stars with sufficient accuracy. The calculation of 
the ^5*0 pairing gaps in pure neutron matter, or symmet- 
ric nuclear matter based on bare NN interaction depends 
strongly on the forces that are used. 

In this work, we compare the results of the zero-range 
and finite-range interactions as pairing interactions. For 
the finite-range interactions, the effective Gogny DIS [l^l 



141 interactions are used to calculate the 



or M3Y-P4 
pairing field. 

For the zero-range force in the pairing channel, we use 
the form 

Vir,~r2) = Vo(^l-v(^^^ ^ir, ^ r^) (16) 

where Vq, tj and a are parameters that can be adjusted. 
We have studied the predictions of two types of pair- 
ing forces: a density-independent interaction (77 = 0) 
that gives rise to volume pairing and a density-dependent 
delta force that can give rise to surface pairing. In Ap- 
pendix B, we show the details of the calculations of pair- 
ing matrix elements. 

Up to now, the magnitude of pairing correlations in 
neutron matter is still a subject of debate. The DIS 
Gogny interaction commonly used in finite nuclei calcu- 
lations gives a maximum value of pairing gap in infinite 
neutron matter of about 2.4 MeV at a Fermi momentum 
kp ~ 0.8 fm^^ (22l |. The maximum value of pairing gap 
in infinite matter using the M3Y-P4 interaction is about 
3.2 MeV at the similar Fermi momentum |14| . On the 
other hand; the microscopic calculations of Refs. (23l.[2^ 
predict for the maximum gap a value of about 1 MeV. 
One observes that there are three different scenarios for 
pairing correlations in neutron matter. In our study, we 
use the density dependent delta force for the pairing in- 
teraction to simulate the third scenario. This is obtained 
by adopting the parameter values of Ref.Q: Vb=-330 
MeV fm-3, 7?=0.7, and a=0.45 .Thus, for each WS cell 



TABLE I: The Wigner-Seitz cells considered in this work, p, 
N, Z and Rws are the baryonic densities, the numbers of 
neutrons and protons, and the WS cell radii, respectively. All 
values are taken from Ref. [l|]. 
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-Rws [fm] 
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0.143 
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10" 
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140 


40 


53.6 
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0.250 
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10" 


2 


160 


40 


49.2 
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0.375 
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10" 


2 


210 


40 


46.3 


7 


0.549 


X 


10" 


2 


280 


40 


44.3 


6 


0.994 


X 


10" 


2 


460 


40 


42.2 


5 


0.233 


X 


10" 


1 


900 


50 


39.3 


4 


0.361 


X 


10" 


1 


1050 


50 


35.7 


3 


0.557 


X 


10" 


1 


1300 


50 


33.1 


2 


1.275 


X 


10" 


1 


1750 


50 


27.6 


1 


2.968 


X 


10" 


1 


1460 


40 


19.6 





4.931 


X 


10" 


1 


950 


32 


14.4 



we perform three HF-BCS calculations with three differ- 
ent pairing forces. The BCS pairing window is chosen to 
be ± 6 MeV around the Fermi level. 



III. RESULTS AND DISCUSSIONS 

In the present study, we have performed the HF-BCS 
calculations for a set of 11 representative WS cells de- 
termined in Ref. [if. The considered density range is 



from neutron drip density pr, 



1.743 X 10"^ po to 



about pmax = 0.5po- In this densit y ra nge, the nuclear 
clusters are assumed to be spherical [13] • Above the den- 
sity Pmax the energy per baryon approaches the value 
of the uniform neutron system and the cells in the in- 
ner crust might deviate from the spherical shape [27| . 
Following Ref. [l| we denote the WS cells like a nu- 
cleus with Z protons and TV neutrons. The eleven zones 
of the representative cells of the inner crust with mean 
densities and corresponding proton number Z and neu- 
tron number N in each cell are listed in Table HI The 
decreasing zone number Nzone are from 10 to 0, corre- 
sponding to the increasing density from the minimum 
density pmin = 2.79 x 10~^ fm~'^ to the maximum den- 
sity Pmax = 7.89 X 10"^ fm~^. The WS cells are denoted 
like a nucleus as i^^Zr, ^oozr, 250zr, 3202^, 5002^, gsogn, 
"OOSn, isooSn, i^sosn, i5002r, and ^s^Qe, as introduced 
in Ref. [l|. These (N,Z) values must be considered as 
indicative since they were determined by using a differ- 
ent model for the HF mean field in Ref.(l|. The WS cell 
radii Rws Eire calculated by the following relation 



<p>-- 



A 



47r p3 
3 ^WS 



(17) 



where p and A are the density and mass number of the 
considered cell, respectively. The values of the radii Rws 
are also shown in the last column of Table HI 

In isolated finite nuclei, for a given number of protons 
there is always a maximum number of bound neutrons. 
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This neutron stability limit defines the neutron drip line. 
Since the neutron-rich nuclei are quickly beta decaying, 
then the neutron drip line is usually drastically limited in 
the laboratory. This is not the case for the neutron-rich 
systems immersed in the inner crust of neutron stars. In 
the case of WS cells, the beta decay is blocked by the 
presence of the degenerate electron gas uniformly dis- 
tributed inside the cell. Therefore, in this case the nuclei 
inside the inner crust of neutron stars can bind more neu- 
trons than the nuclei in the vacuum. Furthermore, there 
are many delocalized neutrons forming a uniform neutron 
gas and filling the outer region of the WS cells. 



A. HF and HF-BCS density distributions 

First, we discuss the case of the WS cells calculated 
in the HF approach without pairing effects. Fig. [T] dis- 
plays the HF proton and neutron density profiles of ^^'^Zi, 
2"0Zr, ^sozr, ^^^Zi and ^°^Zt systems obtained with DIS, 
M3Y-P4 and SLy4 interactions. One notes that the nu- 
merical HF calculations with SLy4 are the same as in Ref. 
. We observe that the HF calculations with the finite- 
range DIS and M3Y-P4 interactions give very similar 
results, while the results obtained with SLy4 interaction 
are different. 

We will analyze the case of the cell ^^°Zr. One knows 
that the depopulation of the s state can lead to a deple- 
tion of the central density [2^, [2^. In the case of this 
cell, the 5si/2 neutron state is fully filled with SLy4 in- 
teraction and empty with DIS and M3Y-P4 interactions. 
Thus, the neutron densities in the center of the cell ob- 
tained with DIS and M3Y-P4 interactions are smaller by 
a factor of 2.3 than that obtained with SLy4 interaction. 
However, the neutron gas density obtained with DIS and 
M3Y-P4 interactions is more than 1.8 times greater than 
that of SLy4 interaction. We have checked that if one in- 
tegrates the neutron density up to r = 10 fm, where the 
neutron density profile becomes approximately constant, 
one finds about 80 neutrons and 90 neutrons with DIS 
(or M3Y-P4) and SLy4 interactions, respectively. Thus, 
it appears that large surface regions are observed in the 
case of DIS and M3Y-P4 interactions, and the neutron 
gas density of the cell ^^"Zr is much higher. Similar sit- 
uations happen in two other cells ^"^Zr and ^^°Zr. In 
the case of cells ^^"Zr and ^"^Zr, where the 5si/2 neutron 
state is fully filled, one observes that the nuclear cluster 
region becomes larger with the SLy4 interaction, while 
its neutron gas density is smaller by a factor of 2 than 
those obtained with DIS and M3Y interactions. For the 
above cells one concludes that the neutron density in the 
center of the cell becomes smaller with DIS and M3Y-P4 
interactions. 

However, the situation is opposite at higher densities, 
such as in the cells ^^°Sn, "O^Sn, i8°°Sn, ^^^°Sn and 
1500 2].^ These changes can be seen in Fig. [2J where the 
neutron densities of the nuclear clusters obtained with 
SLy4 interaction are always smaller than those obtained 



with DIS or M3Y-P4 interactions. One can also see that 
the surface thickness of the nuclear cluster with SLy4 in- 
teraction becomes larger by about 10% in the cells ^^°Sn, 
"OOSn, isoogn and ^^^°Sn. In the cell ^^°°Zr, the nu- 
clear cluster surface is similar with the three interactions. 
Since the neutron density in the center of this cell calcu- 
lated with SLy4 interaction is smaller than that obtained 
with DIS or M3Y-P4 interactions, then its outer neutron 
gas density is a little higher. For the highest density cor- 
responding to the cell ^^^Ge, there is still some trace of a 
central cluster and of an outer neutron gas in the case of 
SLy4 whereas this separation fades away with DIS and 
M3Y-P4 interactions. 

As mentioned above, the first microscopic calculations 
of the WS cells in the inner crust matter were done in 
Ref. The pairing effects were not taken into ac- 

count in these calculations because it was assumed that 
their contributions would be small in comparison with 
the total binding energy of the considered system. How- 
ever, calculations using an energy functional involving 
the neutron and proton pairing correlations by Baldo et 
al. show that the mean density p of the equilibrium 
configuration (Z, Rws) can be changed significantly due 
to the pairing effects. We therefore calculate also the WS 
cells in HF-BCS approximation. The DIS and M3Y-P4 
interactions are used to construct the mean field and the 
pairing field. For comparison we have also considered a 
hybrid case where a density-dependent delta force is cho- 
sen for the pairing interaction and the DIS interaction is 
used in the mean field channel. 

Fig. [3]shows the proton and neutron densities obtained 
in HF and HF-BCS ap pro ximations for 6 WS cells, in 
which the DIS Gogny [ij] interaction is used in both 
the mean field and pairing channels. One can see that 
the behaviour of the proton and neutron densities ob- 
tained in the two approximations are similar for ^^''Sn 
and ^^'^"Sn cells corresponding to the range of density 
p - 2.79 X 10-4 fm-3 to 5.77 x lO'^ fm'^. One can 
conclude that the pairing effects are very small on these 
cells. For the higher density region p ^ 8.91 x 10"'^ 
fm~^ to 4.75 x 10^^ fm~^, corresponding to cells ^^^''Sn, 
isoogj-^ g^jjj i5002i-^ there are differences between the den- 
sity profiles due to the pairing effects. This is because the 
occupancy of the 8si/2 state is modified in HF-BCS cal- 
culations. Thus, the difference between HF and HF-BCS 
neutron densities at r=0 is around 10-11% for the three 
cells above. Indeed, the occupancies of the 8si/2 neutron 
orbital are 0.91, 0.98 and 0.95 for cells "50Sn, i^oog^ .^^^ 
15002]- in HF-BCS, respectively, while this state is fully 
filled in HF calculations of these cells. In spite of the cell 
isooZr, the HF-BCS neutron densities of the ceUs "50Sn, 
isoog]-] i^ave an extended "surface" before they reach the 
constant values corresponding to the neutron gas. The 
neutron gas densities of these cells are similar in both 
HF and HF-BCS approaches. For the highest-density 
cell ^s^Ge with p =7.89 x 10"^ fj^-s^ ^j^g behaviour of 
the neutron density is slightly changed due to the pairing 
effects. Although the DIS interaction is used in the pair- 
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FIG. 1: (Color online) The HF proton and neutron densities in zones 10 to 6 calculated with the Dirichlet-Neuman boundary 
conditions and using DIS [l3[, M3Y-P4 [li| and SLy4 ^ interactions. 
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FIG. 2: (Color online) Same as Fig. [U but for higher density zones (from zone 5 to zone 0). 
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ing channel, it cannot produce a constant density around 
the outer edge of this ceU. It seems that the cell ^®^Ge 
most probably belongs to the deformed pasta phase. 

Since it is not yet well established what are the pair- 
ing properties of neutron matter, we perform again the 
HF-BCS calculations for all WS cells of the inner crust 
matter using three different pairing interactions, namely 
the DIS, M3Y and delta interactions. The results of pro- 
ton and neutron density distributions obtained with these 
three pairing interactions are shown in Fig. 21 Different 
combinations of mean fields and pairing interactions are 
used in this figure. One observes that the main features 
of the WS cells can be obtained with three kinds of pair- 
ing interactions. At the high density region, except for 
the cell ^■^^^Sn where the pairing effects on the neutron 
density are strong, the calculated proton and neutron 
densities are similar for the other cells such as ^^™Sn, 
isooZr, isooSn and ^^^Ge. In the case of the ^^^°Sn ceU, 
the neutron density distribution of the nuclear cluster cal- 
culated with the zero-range pairing force is higher than 
those obtained with DIS or M3Y-P4 pairing interactions 
because the occupancy of the S>si/2 neutron orbital cor- 
responding to the delta-pairing interaction is larger than 
those of the DIS and M3Y-P4 pairing interactions. In all 
the cells we observe that the density distribution of the 
nuclear clusters obtained with the M3Y-P4 interaction is 
slightly less extended than those obtained with the DIS 
interaction. This effect may come from the differences 
of the ranges and the values of pairing gaps in infinite 
matter corresponding to these two effective interactions. 

B. Pairing fields 

Finally, we discuss briefly the BCS neutron pairing 
fields in the case of a delta -pairing interaction. Rewrit- 
ing Eq. (dni as 

V{r,-r2) = Veff{p{r))S{r,-r2) , (18) 

the pairing field A(r) is a local function which can be 
expressed in terms of the local pairing density K{r) as 

A(r) = V,ff{p{r))K{r) , (19) 

where the BCS pairing density K{r) is 

i 

Here, the factors and Vi are the usual BCS amplitudes. 

The results calculated for the WS cells from zone 10 
to zone 1 (using DIS in the particlc-holc mean field) are 
shown in Fig. [5l except for zone which belongs to the 
deformed pasta phase. In the upper panel of Fig. [5l one 
can see that the pairing field is becoming very small al- 
most everywhere in the cell. However, one can observe 
that in passing from the low density region of the neu- 
tron gas towards the higher density region of the cluster, 



the pairing field is increasing in the intermediate den- 
sity region of the cluster surface in the cells '^^°Zr and 
^"^Zr. This is a manifestation of the bell shape depen- 
dence of the pairing gap on density. For the cells corre- 
sponding to high baryonic densities shown in the lower 
panel of Fig. [5] the slope of the pairing field is changing 
very slowly when it is crossing the region between the 
nuclear cluster and the uniform neutron gas. Using the 
same density-dependent delta force in the pairing channel 
and the DIS force in the mean field channel, the pairing 
fields of the cells ^^^Sn, ^^o^Zr and ^^°°Sn are about three 
times smaller (in absolute value) than those obtained by 
the HFB calculation using the SLy4 interaction in the 
mean field (see Figs. 5 and 6 in Ref. 3). We see that 
the pairing correlations are reduced strongly for all the 
cells in the present work, due to the use of finite range 
interactions in the mean field. 

One can conclude that the behaviour of the pairing 
field in the inner crust matter is rather complex. On the 
other hand the magnitude of the pairing field inside the 
inner crust depends strongly on the scenario used for the 
pairingproperties of infinite neutron matter as shown in 
Refs. 11 g. 



IV. CONCLUSIONS AND PERSPECTIVES 

In this work, we have studied the properties of the 
WS cells in the inner crust matter of neutron stars using 
finite-range, density-dependent interactions in HF and 
HF-BCS approximations. Calculations arc performed 
for 11 representative WS cells by imposing Dirichlet- 
Neumann boundary conditions at the edge of the cell. 
The study is done with the DIS and M3Y-P4 effective 
interactions in both mean field and pairing channels. For 
the pairing field we have also used a density-dependent, 
zero-range force whose parameters have been fixed to re- 
produce the pairing properties of infinite neutron matter 
as predicted by microscopic calculations which take into 
account polarization effects For the HF-BCS cal- 

culations using the zero-range pairing force, the HF mean 
field is calculated with the DIS interaction. 

With the three different pairing interactions (DIS, 
M3Y-P4 and delta force), it is found that the behaviour 
of the proton and neutron density distributions are very 
similar in the low density region, p ^ 2.79 x 10"''' fm~'^ 
to p ~ 5.77 X 10"'^ fm^-^. One can conclude that the 
pairing effects are very small in this region. In the higher 
density region, the pairing effects can make an extended 
"surface" for tire neutron density distribution in the cell 
before they reach a constant density around the outer 
edge of the cell. It should be noted that pairing effects 
are weakened by the use of finite range interactions in 
the particle-hole (mean-field) channel, compared to the 
use of Skyrme interactions. However, one cannot ob- 
tain a constant neutron gas density in the outer part of 
the ^^^Ge cell with finite-range density-dependent inter- 
actions. This cell corresponds to the highest density p 
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FIG. 3: (Color online) The proton and neutron densities obtained with the DIS [13 
mations. The calculations are done from zone 5 to zone 



interaction in HF and HF-BCS approxi- 



= 7.89 X 10 ^ fm and it seems not to belong to the 
spherical case as assumed in our study. 

In all the cells we observe that the density distribution 
of the nuclear cluster obtained with the case of M3Y-P4 
interaction is slightly smaller than those obtained with 
the DIS interaction. This effect may come from the dif- 
ference of the ranges and the values of pairing gaps in 
infinite matter given by the two kinds of effective inter- 
actions. The largest effect on the densities comes from 
the use of finite range interactions in the "particle-hole" 
mean field, compared to the use of Skyrme interaction. 
This result shows that finite range interactions should be 
considered in order to describe microscopically Wigncr- 
seitz cells. 



Appendix A: Hartree-Fock potentials with finite 
range interactions 

We give here some general expressions of HF poten- 
tials corresponding to Gogny-type forces. They can be 
easily adapted to the case of the M3Y-P4 Nakada's in- 
teraction jl4| if one adopts a zero-range approximation 
for the two-body spin-orbit component, as mentioned in 
subscc. II. A . 

The contributions of the central force arc separated 
into a Hartree (direct) and a Fock (exchange) potential: 

j u=l 

X / u2(r2)w^(ri,r2)r2dr2 , (Al) 
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Ufiri,r2) = '^'^jj Uj{ri)uj{r2)vl{ri,r2) 

jL v=l 

I J LV\^2^2jk k 1/2 \' 
,0 Oj L ] 

X (W.Sg^g^ - H„) + Bjg^g^ ^ M , (A2) 

with |/i — /jl < L < {li + Ij), and we use the notation 
J = i2j + 1)1/2. 
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FIG. 4: (Color online) The HF-BCS proton and neutron densities obtained in the three cases indicated in the panels. The 
notation "Gogny-Gogny" means the DIS Gogny interaction is used in the mean field and pairing channels, and similar notations 
for other cases. The calculations are done from zone 5 to zone 0. 



The density-dependent component of the force is sim- 
ilar to the zero-range density-dependent term of Skyrme 
interactions [1, , and we recall for completeness the 
corresponding contribution (direct plus exchange) to the 
mean field: 



K^'i^) = |{(2 + ^3)(2 + «)p"+^(r) - (2x3 + 1) 
x[2p"{r)p,{r) + ap"-\r){pl{r)+pl{r))]} , 

(A3) 

where ia, X3 and a are the parameters of the density- 
dependent force, p is the total nucleon density and q 
stands for protons or neutrons. 

The direct plus exchange spin-orbit mean field is Q : 



are: 



J G protons 



Pp{r2)vo{ri,r2)rldr2 



(A5) 



Vf'^iri,r2) = e^Xl^'?..-i/2'^9.,9.^«^^/-?'/"i(^i)"j(^2) 



where w£(ri, r2) — r^/r^^^ are the multipoles of l/|ri 
|r2|. 



V^ir ) 



-W,i^~{p{r)+p,{r))l.s 
\j{r) + j'{r)^\j,{r) + j'^{r) 



= W^'''\r)l.s + W;^'=\r) , 
where J(r) is the spin density, and J' = 



(A4) 



Finally, the direct and exchange Coulomb mean fields 



Appendix B: Pairing matrix elements 

Here, we give the main expressions for calculating the 
pairing fields in the BCS approximation. The gap equa- 
tion is Ei: 



(Bl) 
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FIG. 5: (Color online) The neutron pairing field A(r) of Eq.fTO} in zones 10 to 1, calculated witli the zero range pairing 
interaction and the DIS interaction in the particle-hole mean field 



where 



Goiaabb) aa\Vp{l,2)\bb 

>oo 



(B2) 



is the J = particle-particle matrix element of the 
Vp{l,2) pairing interaction. 



1. Particle-particle matrix elements with a 
finite-range interaction 



For any interaction of the DIS or M3Y-P4 type the 
general particle-particle matrix element of the finite 
range part is: 



< ac\Vp{l,2)\bd>jM 
J rjdnridr, ^ (-1)^ 



uLSK 



X A,{S)vl{ri,r2)Ra{ri)Rb{ri)R,{r2)Rd{r2) 
' 3d jb J 

. ja jc K 



(B3) 



Here, the total phase is T = L + 5 + /\ + J + ^5 + 
Ay{S) is given by Eq. Ri{r) is the radial part of the 
single-particle wave function of Eq. and 



< > is the reduced matrix element: 



z S L'Jlilj 




(B4) 



For the pairing matrix elements < aa\Vp{l,2)\bb >oo 



we thus obtain: 



< 



aa|yp(l,2)|66 >oo = Ja ^ k ^ i rjdnrldr^ ^ A^{S) 



X V 



l{rur2)Ra{ri)Rb{ri)Ra{r2)Rb{r2) 
X (-1)^+^1 <ya||T(^^^^||y6> |'(B5) 

The zero-range part of the Gogny or M3Y-P4 interac- 
tions does not contribute to the pairing matrix elements, 
while the small contribution of the zero-range spin-orbit 
component to the pairing field is neglected. 



2. Particle-particle matrix elements with a 
zero-range interaction 

The density-dependent delta interaction is taken of the 
form: 



PQ 



<5(ri - ra) (B6) 
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The delta function is expanded on the basis of spherical where 
harmonics: 



<5(ri - r^) = Y.^-ir^^{l)Y-^{2) (B7) 

laab, = Vo [ (^) ) Rl{T)Rl{r)dr . (B9) 

The calculation of the coupled matrix element is similar V \ Po / / 



to that of subsection B.l. The result is: 



< aa\Vp{l,2)\bb >oo= ^(-1)'"+'" ^/aahb , (B8) 
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